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Complete  elliptic  integral  of  the  second  kind 
Distance  factor 

Complete  elliptic  integral  of  the  first  kind 

Modulus  of  elliptic  function,  Equation  (9) 

Radius  vector 

Unit  vectors 

Propeller  radius 

Radius  of  ring  vortex 

Cylindrical  co-ordinates 

Vector  tangential  to  the  vortex  ring 

Induced  velocity 

Axial  induced  velocity 

Rectangular  coordinates 

Circulation 

Advance  ratio 


ABSTRACT 


NTh£s  report  gives  the  axial  induced  velocity  ahead  of 
an  infinitely  bladed  propeller.  The  propeller  is  simulated 
by  a  close  succession  of  ring  vortices  whose  strength  vary 
along  the  propeller  radius  and  which  extend  from  the  propel¬ 
ler  plane  to  Infinity.  The  results  are  compared  with  those 
obtained  for  a  propeller  represented  by  a  uniform  sink  disc. 

INTRODUCTION 

The  flow  field  in  the  vicinity  of  a  propeller  can  be 
calculated  by  replacing  the  propeller  by  means  of  flow  singu¬ 
larities.  The  usual  representation  has  been  to  simulate  an 
infinitely  bladed  propeller  in  a  steady  inviscid  incompressible 
fluid  by  means  of  a  distributed  sink  distribution  over  the 
entire  propeller  disc.  The  problem  was  originally  formulated 
by  Dickmanl  for  a  propeller  having  a  uniform  distribution  of 
loading  along  the  radius  and  the  induced  velocities  resulting 

p 

for  such  a  distribution  has  been  given  in  the  literature  » 
Furthermore,  a  propeller  with  a  radially  varying  distribution 
has  also  been  simulated  by  superimposing  sink  discs  with  a 
common  center^. 
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The  approach  of  the  prerent  method  has  been  to  simulate 
an  infinitely  bladed  propeller  with  radial,  y  varying  distri¬ 
bution  of  loading,  by  means  of  a  series  of  ring  vortices  ex¬ 
tending  from  the  propeller  to  infinity.  Such  a  representation 
1  not  expected  to  give  radically  different  results  from  the 
sinX  representation  for  relatively  large  distances  from  the 
propeller  plane.  In  the  vicinity  of  the  propeller  plane, 
however*  the  magnitude  of  the  induced  velocities  will  be  sub¬ 
stantially  affected  by  the  mathematical  representation  of  the 
propeller. 

The  present  report  is  confined  to  the  formulation  and 
calculation  of  the  Induced  axial  velocities  for  a  propeller 
having  an  optimum  distribution  of  loading.  Similar  formula¬ 
tions,  however,  can  also  be  performed  for  both  the  tangential 
and  radial  induced  velocities. 

Analysis 

The  induced  velocities  in  the  vicinity  of  a  propeller 
are  a  result  of  the  system  of  free  vortices  which  originate  at 
the  propeller  blades  and  constitute  the  propeller  slipstream. 

A  propeller  with  a  very  large  number  of  blades  (i.e.,  an 
infinite  number  of  blades)  can  be  considered  as  shedding  a 
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system  of  closely  spaced  helical  vortex  surfaces  which  are 
equivalent  to  a  compact  system  of  helical  vortex  sheets.  These 
helical  surfaces  can  be  resolved  into  two  parts;  (1)  a  close 
succession  of  transverse  vortex  rings  with  center  at  the  propel¬ 
ler  axis  and  (2)  a  system  of  vortex  lines  which  are  parallel  to 
the  propeller  axis  (Figures  la  and  lb).  The  axial  induced 
velocity  at  any  point  in  space  is  a  result  of  only  the  vortex 
rings ,  as  the  parallel  vortex  lines  do  not  induce  an  axial 
component.  Similarly,  the  tangential  induced  velocities  are 
induced  only  by  the  parallel  vortex  lines. 

In  order,  therefore,  to  evaluate  the  axial  velocities  we 
will  consider  a  vortex  ring  of  radius  r’  with  the  coordinate 
system  r,#-,  z  and  x,y,z  with  z  in  the  direction  of  the  axis 
of  the  ring  as  shown  in  Figure  2. 


By  the  law  of  Biot-Savart 


dv 
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giving  the  velocity  vector  dv  induced  by  an  element  of  length 
ds  -  r'd©  *  at  a  point  (r,0,z)  or  (x,y,z).  With  the  unit 
vectors  T,3\W  in  the  direction  of  x,y,z  axes  we  can  write 

.d^L  —  -  1  +  j  +  °  ■  ^ 
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and  F  is-  the  radius  vector  from  the  vortex  element  to  the 
point  P. 


p  =  ttoiQ'y*.  t  Jp  +  2  k 

and 

P  =  /  z\  l 

(2) 

Hence , 

Px  ds  »  2c«s  B1  i  *■  z  $'>*&'  J  -  [  iCoi(6-9')  -  fe.'1  jc 

The  z-component  of  the  induced  velocity  will  be  obtained  by 
integrating  the  entire  ring. 

yJv).-J_  \  -j- —  (9~ -9  -.•■=z=Tt3  de' 

j  y  zS  V*-  2.tVc*s(*-e0j 

(3) 


This  integration  can  be  performed  by  means  of  complete 
elliptic  integrals  (Ref.  4)  to  give 
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where 


(5) 


and  where  K(k)  and  E(k)  are  the  complete  elliptic  integrals 
of  the  first  and  second  kind. 


For  a  vortex  ring  placed  at  the  axial  position  z0, 
Equation  (4)  becomes 


where 


For  a  propeller  having  an  infinite  number  of  blades,  the 
optimum  circulation  distribution  along  the  radius  is  given  by 

r  =  c 

Where  C  is  a  constant,  As  is  the  advance  ratio  and  E  is  the 
propeller  radius. 

In  order,  therefore,  to  find  the  induced  axial  velocity 
for  the  entire  system  of  ring  vortices,  Equation  (6)  can  be 
written  as: 
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where 


and 


i  =  r'/E 


d 


Z0/H 

z/E 

4  it*/*) 


(9) 


The  induced  velocity  at  any  point  (r/E,d)  is  given  by 
Eauations  (8)  and  (9).  This  velocity  can  be  normalized  in 
terms  of  the  axial  induced  velocity  at  the  propeller  disc. 

For  points  ahead  of  the  propeller,  we  can  write 

1 

“*■  *  ‘  ’Va  l*/*.,  0)  (io) 

where  ga  is  known  as  a  distance  factor  for  the  axial  velocity. 
Similar  expressions  can  easily  be  obtained  for  the  tangential 
induced  velocity  by  considering  the  system  of  bound  vortices 
and  line  vortices. 
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Table  I 


r/: 

d\ 

1.0 

2.0 

3,0 

0.02 

0.032 

0.036 

0.039 

0.10 

0.156 

0.172 

0.186 

0.20 

0.3° 

0.295 

0.412 

IM 

0.339 

0.462 

0.40 

0.50 

0.509 

0.588 

0.537 

O.615 

0.559 

0.635 

0.60 

0.652 

0.678 

O.696 

0.70 

0.705 

0.728 

0.744 

O.9O 

0.783 

0.801 

0.818 

1.10 

0.833 

0.851 

0.861 

1.50 

0.899 

O.9O9 

0.915 

1.90 

0.933 

0.939 

0.944 

=  0.0 

4.0 

0.041 

0.197 

0.355 

0.478 

0.574 

0.649 

0.708 

0.755 

0.822 

0.867 

0.919 

0.946 


5*0  6,0 


0.044 

0.2064 

0.367 

0.4Q0 

0.585 

O.659 

0.717 

0.763 

0.829 

0.872 

O.922 

0.948 


0.047 

0.214 

0.377 

0.500 

0.594 

0.667 

0.724 

0.769 

0.833 

0.875 

0.Q24 

O.950 


0.02 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

O.9O 

1.10 

1.50 

1.90 


0.036 

0.168 

0.309 

0.425 

0.521 

0.598 

0.660 

0.711 

0.787 

0.838 

O.899 

0.933 


0.040 

0.184 

0.331 

0.452 


0. 


0.622 

O.683 

0.731 

0.803 

0.851 

0.908 

0.939 


r/K  =  0.2 


0.043 

O.I96 

o°:$ 

0.564 

O.638 

0.697 

0.745 

O.813 

0.859 

O.913 

0.942 


0.045 

0.203 

O.360 

0.480 

O.574 

0.648 

0,706 

0.752 

O.8I9 

0.864 

°-9l7 

0.944 


0.047 

0.209 

O.366 

0.487 

0.581 

0.654 


0.712 

0.757 

0.823 

O.867 

O.919 
O'.  946 


0.047 

0.212 

O.37I 

0.492 

0.586 

0.658 

0.715 

0.761 

0.826 

O.869 

0.920 

0.947 
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Table  I  (Cont.) 
r/R  =0.3 

2.0  3.0  4.0 


5.0  6.0 


j 


0.02 

O.O38 

0.041 

0.043 

0.10 

0.176 

0.190 

0.197 

0.20 

O.3O 

O.32O 

0.436 

0.339 

0.457 

S:& 

0.40 

0.529 

0.550 

0.5bl 

0.50 

0.605 

0.624 

0.634 

0.60 

0.665 

O.683 

0.693 

0.70 

0.715 

0.731 

0.739 

O.9O 

0.788 

O.801 

O.808 

1.10 

0.838 

0.849 

0.854 

1.50 

O.899 

O.9O6 

0.910 

1.90 

0.932 

0.937 

0.939 

0.044 


0.200 

0.354 

0.473 

0.566 


0.639 

O.697 

0.744 

0.812 


0.857 

0.912 

0.941 


0.045 

0.201 


0.356 

0.476 


0.569 

0.642 


0.700 

0.746 

0.814 


0.859 

0.913 

0.941 


0.044 

0.201 

0.357 

0.477 

0.571 

0.643 

0.702 

0.748 

0.815 

0.860 


r/R  =  0.4 


0.02  0.040 
0.10  0.184 

0.20  0.330 
O.3O  0.445 
0.40  0.537 
0.50  0.610 
0.60  0.668 
0.70  0.716 
0.90  0,787 
1.10  0.836 
1.50  O.896 
1,90  O.930 


0.042 

0.043 

0.192 

0.194 

0.341 

O.343 

0.457 

0.460 

0.548 

0.552 

0.621 

0.624 

0.679 

0.682 

0.726 

0.729 

8:81 

0.798 

0.845 

0.901 

0.933 

0.903 

0.934 

0.042  0.041 
0.193  0.191 
O.343  0.342 
0.460  0.459 
0.553  0.552 
0.625  0.625 
O.683  0.683 
0.730  0.730 


0.903  0.904 

0.935  O.935 


0.042 

0.191 

0.341 

0.459 

0.552 

0.625 

O.683 

O.73O 

0.800 

0.847 

0.904 

0.935 
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Table  I  (Cont.) 


r/K  =0.5 


2.0 


3.0  4.0 


5.0 


0.02 

0.10 

0.20 

O.3O 

0.40 

0.50 

0.60 

0.70 

O.9O 

1.10 

1.50 

1.90 


0.042 

0.191 

0.337 

0.451 

0.539 

O.6O9 

0.662 

0.712 

0.782 

0.830 


O.89I 

0.925 


0.042 

0.041 

0.191 

0.188 

0.337 

O.334 

0.451 

0.447 

0.540 

0.537 

0.611 

0.608 

0.66  7 

0.665 

0.714 

0.712 

0.783 

0.782 

0.832 

0.831 

O.892 

0.892 

0.926 

O.926 

0.040 

0.185 

°*33° 

0.443 

0.533 

0.605 

O.663 

0.710 

0.781 

0.830 

O.891 

0.926 


0.040 

O.I83 

0.327 

0.441 

0.531 

O.6O3 

0.661 

0.708 

0.780 

0.82Q 

O.89I 

0.925 


6.0 

O.O39 

0.181 

0.325 

0.438 

0.529 

0.601 

0.659 

O.707 

0.77? 

0.828 

O.89O 

O.925 


\ 

.1 


i'/K  =  0.6 


0.02  0.043 

0.10  0.194 

«  A  A  ' 


0.40  .O.532 
0.50  0.600 
0.60  0.654 
0.70  O.699 
0.90  0.768 
1.10  0.817 
1.50  0.880 
1.90  O.9I7 


0.041 
0.187 
0.328 
0.436 
0.521 
O.589 
0.645 
0.691 
0.762 
0.812 
0.8  77 
O.915 


0.040 
0. 
o._ 
0.425 
0.511 
0.580 
0.637 
0.684 
0.756 
0.808 
0.874 

O.913 


0.039 

0.175 

O.311 

0.418 

0.504 

0.574 

O.631 

O.679 

0.752 

0.805 

0.873 

O.912 


0.037 

O.172 

0.307 

0.414 

0.500 

0.570 

0.628 

0.675 

0.750 

0.803 

0.871 

0.911 


0.037 

0.170 

0.304 

0.411 


OA97 

0.567 

0.625 

0.673 

0.748 

0.802 

0.870 

O.9IO 
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Table  I  (Cont.) 


r/R  =0.7 

2.0  3.0  4.0 


5.0  6.0 


0.02  0.043 
0.10  0,191 
0.20  0.326 
0.30  0.427 
0.40  0.506 
0.50  0.570 
0.60  0.623 
0.70  0.668 
0.90  0.738 
1.10  0.790 
1.50  0.859 
1.90  0.901 


0.174  0.163 
0.302  0.286 
0.401  0.384 
0.480  0.463 
0.546  0.530 
0.601  0.586 
0.648  O.634 
0.722  0.711 
0.777  0.768 
0.851  0.845 
0.8§5  O.89I 


0.034  0. 

0.157  0. 

0.277  0. 

0.373  0. 

0.453  0. 

0.520  0.514 

0.577  0.571 

0.626  0.621 

0.704  0.700 

0.763  0.760 

0.841  0.839 

0.888  O.887 


0.151 

0.268 

o°$I 
0.510 
0.56  7 
O.617 
0.697 
0.757 
0.837 
0.886 


r/R  =  0.8 


0.02 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.90 

1.10 

1.50 

1.90 


0.039 

0.162 

0.274 

0.358 

0.429 

0.490 

0.544 

0.592 

0.670 

0.730 

0.814 

0.866 


0.032 

0.135 

0.232 

0.312 

O.382 

0.446 

0.502 

0.553 

0.637 

0.703 

0.795 

0.853 


0.028 

0.120 

0.209 

0.285 

0.355 

0.419 

0.477 
0.530 
0.618 
0.68  7 
0.784 
0.845 


0.026 

0.111 

O.I96 

0.270 

0.340 

"0.404 

0.463 

0.517 

0.607 

O.678 

8:8! 


0.024  0..  024 
0.106  0.103 
0.188  O.I83 
0.261  0.256 
0. 

0. 

0. 

0.508  0.503 
0.600  O.596 
0.6 72  0.669 
0.774  0.772 
0.838  0.836 


331 

395 

455 


0.325 
0 *3,89 
0.449 
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Table  I  (Cont. ) 


r./R  =  0.9 

2.0  3.0  4.0 


5.0  6.0 


0.02 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

O.9O 

1.10 

1.50 

1.90 


0.006 

0.011 

0.032 

0.080 

0.142 

0.-208 

0.272 

o°:$ 

0.536 

0.671 

0.753 


-0.016 

-0.051 

-0.057 

-0.020 

0.039 

0.107 

0.176 

0.244 

0.367 

0.470 

0.624 

0.726 


-0.022 

-0.079 

-O.O99 

-O.O69 

-0.012 

0.055 

0.126 

0.197 

0.326 

0.436 

0.599 

0.707 


-0.026 

-0.093 

-0.121 

-0.095 

-0.041 

0.026 

O.O99 

0.170 

0.304 

0.416 

0,585 

O.697 


-0.028 

-0.101 

-0.133 

-0.111 

-0.057 

0.010 

0.083 

0.155 

O.29O 

0.405 

0.577 


-0.029 

-0.105 

-0.141 


-0.120 

-0.068 

0.000 

0.072 

0.145 

0.281 


0.397 

0.572 

0.687 


Results 


The  distance  factor  ga,  has  been  computed,  using 
Eouatlons  (8),  (9),  and  (10),  on  the  IBM  704  Computer  for 
values  of  r/K  from  0  to  0.9,  d  from  0  to  1.9  and  from 

1.0  to  6.0.  The  results  are  given  in  Table  I,  computed  to 
an  accuracy  of  0.001. 

The  results  can  be  compared  with  those  obtained  from 
a  sink  disc  representation  of  a  propeller  of  uniform  loading. 
The  comparison  is  made  in  Table  II,  and  It  is  noticed  that  the 
difference  is  considerable,  particularly  in  the  vicinity  of 
the  disc.  This  is  primarily  due  to  relatively  large  change  of 
the  axial  velocity  component  near  the  propeller  plane,  indi¬ 
cating  that  the  approximation  of  uniform  loading  along  the 

Table  II 

Comparison  of  ga  with  uniform  sink  disc 


r/R  =  0.6 

r/R  =0.8 

d 

% = l-° 

=  6.0 

Sink  disc 

'/%  - 1-° 

XL  -  6-° 

Sink  disc 

0.3 

0.443 

0.411 

0.373 

0.358 

0.256 

0.477 

0.5 

0.600 

0.567 

0.536 

0.490 

0.389 

0.617 

0.7 

0,699 

0.673 

0.646 

0.592 

0.503 

0.700 

* 


propeller  radius,  leads  to  inaccuracies  of  the  induced 
velocity  field.  Examination  of  the  results,  shows  that  the 
distance  factor  ga  is  larger  than  that  obtained  from  a  uni¬ 
form  disc  representation  for  propeller  radii  smaller  than 
0.6,  however,  for  propeller  radii  of  0.8  or  larger  the 
distancte  factor  is  smaller.  This  would  indicate  that  although 
the  sink  distribution  representation  may  result  in  the  correct 
total  effect,  the  radial  distribution  of  such  velocity  is  a 
sensitive  function  of  the  radial  distribution  of  loading. 

Conclusions 

Comparison  of  the  axial  induced  velocity  in  the  vicinity 
of  a  propeller  indicates  that  the  loading  distribution  is 
important  for  points  near  the  propeller  plane.  The  velocities 
have  been  calculated  for  an  optimum  loading  distribution  which 
in  general  is  not  too  different  from  th'c*  radial  distribution 
usually  used.  In  specific  instances,  however,  it  may  be 
necessary  to  use  the  actual  circulation  distribution  and 
this  can  most  easily  be  performed  by  calculating  the  induced 
velocity  resulting  from  the  differences  in  loading  between 
the  actual  and  optimum  distribution. 
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